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ABSTRACT 

The Kepler space telescope has opened new vistas in exoplanet discovery space by revealing popula- 
tions of Earth-sized planets that provide a new context for understanding planet formation. Approx- 
imately 70% of all stars in the Galaxy belong to the diminutive M dwarf class, several thousand of 
which lie within Kepler's field of view, and a large number of these targets show planet transit signals. 
The Kepler M dwarf sample has a characteristic mass of 0.5 M & representing a stellar population 
twice as common as Sun-like stars. Kepler-32 is a typical star in this sample that presents us with a 
rare opportunity: five planets transit this star giving us an expansive view of its architecture. All five 
planets of this compact system orbit their host star within a distance one third the size of Mercury's 
orbit with the innermost planet positioned a mere 4.3 stellar radii from the stellar photosphere. New 
observations limit possible false positive scenarios allowing us to validate the entire Kepler-32 system 
making it the richest known system of transiting planets around an M dwarf. Based on considerations 
of the stellar dust sublimation radius, a minimum mass protoplanetary nebula, and the near period 
commcnsurability of three adjacent planets, we propose that the Kepler-32 planets formed at larger 
orbital radii and migrated inward to their present locations. The volatile content inferred for the 
Kepler-32 planets and order of magnitude estimates for the disk migration rates suggest these planets 
may have formed beyond the snow line and migrated in the presence of a gaseous disk. If true, this 
would place an upper limit on their formation time of ~ lOMyr. The Kepler-32 planets are represen- 
tative of the full ensemble of planet candidates orbiting the Kepler M dwarfs for which we calculate an 
occurrence rate of 1.0 ± 0.1 planet per star. The formation of the Kepler-32 planets therefore offers a 
plausible blueprint for the formation of one of the largest known populations of planets in our Galaxy. 
Subject headings: planetary systems — methods: statistical — planets and satellites: formation 
- planets and satellites: detection — stars: individual (KID 9787239/KOI- 
952/Kepler-32) 



1. BACKGROUND 

Before the discovery of exoplanets around main se- 
quence stars two decades ago, models of planet formation 
were based on a solitary example: our own Solar System. 
Despite the discovery of hundreds, if not thousands of 
additional planets in the years since then, observational 
and theoretical efforts have focused on the formation of 
planets around Sun-like stars. However, the Sun is not 
a typical star. Seventy percent of stars in the Galaxy 
are dwar fs of the M spectral cl ass ("M dwarfs" or "red 
dwarfs," Bochanski et al.l 120101 ) . with masses that are 
only ~ 10 — 50% the mass of the Sun, much cooler tem- 
peratures, and different evolutionary histories. These dif- 
ferences likely result in different formation and evolution- 
ary histories for their planets. For example, both Doppler 
and transit surveys have revealed a paucity of gas gi- 
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ant planets around M dwarfs, and a relative overabun- 
dance of planets wit h masses less than that of Neptune 
((Howard et al.|[20r2l) . This is in contrast to the high gas 
giant occurrence rate aroun d stars more massive than the 
Sun (j Johnson et alJl2010al ). These correlations between 
stellar mass and gas giant occurrence are likely a conse- 
quence of the lower disk masses around M dwarfs, which 
result in less raw material available for planet building 
(|Laughlin et alj|2004t Kennedy fc Kenvonll2008t ). 

Roughly 5500 of the 160,000 stars targeted by NASA's 
Kepler mission are M dwarfs with a mass distribution 
skewed toward the high mass end of the spectral class 
(see S I5.1I for details). Of these stars, 66 show at least 
one periodic planetary transit signal. Aside from a 
single outlier— the h ot Jupiter system around KOI-254 
(| Johnson et al.l 120121 ) — the ensemble of 100 planet can- 
didates around these M dwarfs have radii ranging from ~ 
0.5 to 3 i?^ and semima.jo r axes within about a fe w tenths 
of an astronomical unit ([Muirhead et al.ll2012aft . These 
compact pl anetary configurations around the lowest mass 
stars (e.g., iMuirhead et al.ll2012bD offer a number of ad- 
vantages: the signal-to-noise ratio of the transit signals of 
close-in planets is boosted due to the increased nu mber of 
transits per observing period (|Gould et al.ll2003h : transit 
depths are larger for a given pl anet radius allowing detec- 
tions of ever smaller planets (|Nutzman fe Charbonneaul 
2008) ; and the reduced temperatures and luminosities of 
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their host stars lead to equilibrium temperatures compa- 
rable to the Earth's, despite their extreme proxim ity to 
their stars (|Kasting et al.lll993l iTarter et al.ll2007l ). 

The Kepler M dwarf planets have been calculated to 
have an occurrence rate a factor of ~ 3 higher than 
for solar-type stars with occurrence r ates increasing as 
planet size and stell ar mass decrease (jMann "eTaT][20T2l: 
iHoward et al.l I2012D . Microlensing studies also sug- 
gest a high planet occurrence around low mass stars 
(jCassan et al.l[2012T ). However, the large errors on these 
results and the fact that microlensing surveys probe plan- 
ets at larger separations from their host stars than tran- 
sit surveys complicate the comparison between samples. 
The high frequency of small planets around low-mass 
stars is compounded by the fact that l ower mass star s 
are more common than solar- type stars (|Chabrierl l2003). 
Therefore the mechanisms by which small planets form 
around the lowest mass stars determine the characteris- 
tics of the majority of planets currently known to exist 
in our Galaxy. 

Kepler-32 is an M1V star with half the mass and ra- 
dius of the Sun, roughly two-th irds the Sun's tempera - 
ture, and 5% of its luminosity (|Muirhead et al.l l2012af ) . 
While showing 5 distinct transit signals, to date only two 
of the planets have been validated, Kepler-32 b and c, 
from the timing variat ions of their transit signals (TTVs; 
IFabrvckv eraLll2012al ). The Kepler-32 host star and its 
planets are typical of the full Kepler M dwarf sample 
except for the chance alignment of this system with re- 
spect to our line of sight that offers a unique and ex- 
pansive view of its dynamical architecture. To capital- 
ize on this chance alignment we bring to bear a suite of 
ground-based observations using the W. M. Keck Obser- 
vatory an d the Robo-AO system on the Palomar 60-inch 
telescope (B aranec et al.l l2012f ) to validate and charac- 
terize the system in detail. We follow this with an in- 
depth analysis that allows us to place tight constraints 
on the formation and evolution of the Kepler-32 plan- 
ets, which in turn has implications for the formation of 
planets around early M dwarfs in general. 

Since all planet parameters are derived directly from 
the physical characteristics of the host star, we start by 
refining the stellar properties in §[5] High spatial reso- 
lution images and optical spectra are then used by our 
false positive statistical analyses to validate the remain- 
ing three transit signals from the Kepler data in §[3l 
These above analyses provide the foundation for an ac- 
curate characterization of the Kepler-32 planets followed 
by an investigation into the formation and evolutionary 
history of this system implied by its observed architec- 
ture in §|4j We then turn to a discussion of the ensemble 
of Kepler M dwarf planets in §[5] and argue that the for- 
mation pathway deduced for Kepler-32 offers a plausible 
blueprint for the formation of planets around the small- 
est stars, and thus the majority of all stars throughout 
the Galaxy. 



2. Kepler-32 STELLAR PROPERTIES 

The observed stellar properties for Kepler-32 are sum- 
marized in Table [TJ It was originally listed in the Ke- 
pler Input Catalog (KIC) with an effective tempera- 
ture of 3911K, below the 4500 K threshold where the 



TABLE 1 

Observed Stellar Properties 



Parameter 


Value 


Ref. 


a (J2000) 


19:51:22.18 


1 


A { T9nnn\ 




i 


fi a (mas/yr) . . 


-8 


1 


Us (mas/yr) . . 


20 


1 


K P 


15.801 


1 


A v 


0.154 


1 


9 


17.251 


1 


r 


15.913 


1 


V 


16.452* 


2 


J 


13.616 ±0.023 


3 


H 


12.901 ± 0.024 


3 


K s 


12.757 ±0.024 


3 


v rad (kms -1 ) 


-32.5 ± 0.5 


1 
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Boruck l et al.l 
I (120051) ; (3) 



IBrown et all H201U ): f4) this work 
Converted from g and r 



photometric method f or determining ste llar parame- 
ters is deemed reliable (|Brown et al.ll20ill) . Follow-up, 
medium-resolut i on infr ared spectroscopy presented by 
iMuirhead et al.l (j2012aj ) revised and refined the KIC val- 
ues to T eff = 3727t^ 2 and [Fe/H]= 0.04 ± 0.13. 

In this work we supplement these measurements with 
an independent estimation of the stellar mass, radius 
and metallici ty using the bro adband photometric method 
presented by Uohnson et al] ([2012D . This method evalu- 
ates the posterior probability distribution of each stellar 
parameter conditioned on the observed apparent magni- 
tudes and colors of the star listed in the KIC using model 
relationships such as the Delf osse et al.l (200C) mass- 
luminosity relation and the iWest et al.l ( 2005 ) color- 
spectral-type relation. Our Markov Chain Monte Carlo 
analysis yields M* = 0.57 ± 0.06 M©, i?* = 0.53 ± 
0.04 R Q , and [Fe/H]= -0.05 ±0.17. 

We c ombine the results of the above analy s es wit h 
thos e oflBorucki et all (l2011al) , IMuirhead et al.l (|2012al ) , 
and IFabrvckv et al.l (|2012a|) by weighted mean to ob- 
tain the final values used in this work presented in 
Table p We adopt 10% errors on the iBorucki et al.l 
(|2011al) values for stellar mass and radius. For the ef- 
fective temperature, we construct probability distribu- 
tion functions and compute the mean and 1 a errors nu- 
merically to ac c ount fo r the asymmetric errors in the 
IMuirhead et al.l (|201 2a) value. Using the iDelfosse et al.l 
( 2000) K-b&nd mass- luminosity relationship we find that 
our final stellar mass and its uncertainty give a distance 
modulus of 7.407 ± 0.098, corresponding to a distance of 
d = 303 ± 14 pc. 

2.1. Stellar Age 

It is difficult to measure the age of field M dwarfs unas- 
sociated with known moving groups or clusters because 
their global stellar properties change little over the course 
of their main sequence lifetime. However, Kepler-32 has 
a slow rotation rate of P = 37.8 ± 1.2 days measured 
from the modulation of the light curve by star spots 
(jFabrvckv et al.|[2012af) implying that it is old. To esti- 
mate the age of Kepler-32 we use the age equation from 
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IBarnesl((2OT0l) . 

The dimensionless constants kc — 0.646 days Myr -1 and 
kj = 452 Myr day^ 1 are determined observationally, the 
convection turnover time is estimated as r = kj /2-f 2 (B— 
V), where w e use f(B - V) = 0.77(B - V - 0.4) 6 
(lBarnes| | 2007f) and (B — V) = 1.5 (converted from g — r 
using Uester et all 12005ft to give r = 155 days. Using 
the median value of the initial spin period, Pq = 2.81, 
as required to produce the observed rot ation rates for 
» 0.6 Mff stars in the Praesepe cluster (|Agueros et al.1 
120111 ) we derive an age of 2.7 Gyr for Kepler-32. The 
large scatter in observed rotation rates for Praesepe 
members used here to calibrate Pq suggests the age of 
Kepler-32 could be as young as 2.3 Gyr or as old as 
3.7 Gyr. For this work we will assume that the age of 
Kepler-32 is greater than 2 Gyr. 



3. VALIDATION OF TRANSIT SIGNALS 

Two of the Kepler-32 planets (b and c) have been val- 
idated previous to this study u sing signatures of thei r 
mutual dynamical interactions (Fabrvck v et al.l l2012aj) . 
The other three transit signals have hitherto retained 
the status of planet candidates. The probabilities for 
Kepler planet candidates to be astrophysical false posi- 
tives, such as blended stellar eclipsin g binaries, are gen- 
erally low ([Morton fc Johnson! 120 111 ), an d may even be 
negli gible in multiply transiting systems (jLissauer et al.1 
12012ft . Diffraction limited images in the Kepler band 
made with th e Robo-AO system on the Palomar 60- 
inch telescope (Barancc ct al. 2012) show no evidence of 
blended companions for a 5 a contrast of Az « 3.5 at 
0"5 and Az sa 4.5 at 1". High-resolution optical spec- 
troscopy with the HIRES spectrograph on Keck I and 
high-resolution adaptive optics imaging in the near in- 
frared using the NIRC2 camera on Keck II place strin- 
gent constrains on the probabilities of various false pos- 
itive scenarios, allowing a probabilist ic valida t ion of the 
transit signals using the procedure of lMorto"nl (|2012h . 

3.1. Keck I/HIRES Spectroscopy 

We obtained spectra of Kepler-32 at Keck Observa- 
tory usin g the Hig h-Resolution Echelle Spectrometer 
(HIRES, iVogt et al.lll994ft on Keck I with the standard 
observing setup used by the California Planet Survey 
(| Johnson et al.ll2010bl ) which covers a wavelength range 
from 3640A to 7820A. Because of the star's faint visual 
magnitude we used the C2 decker corresponding to a pro- 
jected size of 14"0 x 0"851 to allow sky subtraction and 
a resolving power of R = 55000. We obtained 2 observa- 
tions of Kepler-32 (UT 18 June 2011 and UT 13 March 
2012), both with an exposure time of 700 seconds, each 
resulting in a sig nal-to-noise ratio (SNR) of 20 at 6500A. 

To measure the systemic radial velocity of Kepler-32 
and to check for evidence of a second set of stellar lines 
in the spectrum, we performed a cross-correlation anal- 
ysis with respect to a HIRES spectrum of a similar star 
(HIP 86961: V rad = -28.9 ± 0.4 km s" 1 : iNidever et al. 
2002) using the methodology described by I Johnson et al.l 



(|2010al) . Our analysis yields a systemic radial velocity 
of —32.5 ± 0.5 km s -1 . No second set of lines are evi- 
dent above the noise floor in the cross-correlation func- 
tion, allowing us to rule out stellar companions with 
V > lOkms -1 and U-band brightnesses within 2 magni- 
tudes. 

3.2. Keck II/NIRC2 Adaptive Optics Imaging 

We performed adaptive optics imaging using the 
NIRC2 instrument on the Keck II telescope on UT 24 
June 2011 to rule out false positive scenarios involving 
blended sources. With a Kepler bandpass magnitude 
Kp = 15.80, Kepler-32 is relatively faint for natural 
guide star observations. Nevertheless, we were able to 
close the adaptive optics (AO) system control loops on 
the star with a frame rate of 41 Hz. With sufficient counts 
in each wavefront sensor subaperture, a stable lock was 
maintained for the duration of the observations. 
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Fig. 1. — J and K' band contrast curves derived from our Keck 
II/NIRC2 adaptive optics imaging constrain false positive scenarios 
involving blended sources. 

We acquired two sets of dithered images using the 
NIRC2 medium camera (plate scale = 20maspix~ 1 ) a 
sequence of 9 images in the J-filter (central wavelength 
= 1.25 /im), and a sequence of 9 images in the if -filter 
(central wavelength = 2.12/xm). Each image had an ex- 
posure time of 20 s, resulting in a total on-source inte- 
gration time of 180 s per filter. To process the data, hot 
pixels were removed, the sky-background was subtracted, 
and the images were aligned then coadded. We measured 
the contrast achieved by NIRC2 adaptive optics imaging 
by comparing the peak intensity of the star in the fi- 
nal processed image to the intensity of residual scattered 
starlight at small angular separations. Specifically, we 
calculate the standard deviation of flux, a, within a box 
of size 3x3 FWHM, where FWHM is the point-spread 
function full-width at half-maximum spatial scale (also 
the size of a speckle). The standard deviation is eval- 
uated at numerous locations close to the star and the 
results are azimuthally averaged to estimate the local 
radial contrast profile. Figure [TJ shows the 5 a contrast 
curves based on each reduced AO image. The full in- 
strument field of view is 20" x 20", which corresponds 
to 5 Kepler pixels on a side. No obvious contaminants 
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TABLE 2 
Derived Stellar. Properties 



Parameter KIC 1 M12 2 F12 3 this work Adopted Values 

T off (K) 3911 ±200 3727±5g 2 3900 ± 200 ■•• 37931?° 

M*(M©) 0.49 ±0.05 0.52 ±0.04 0.58 ±0.05 0.57 ±0.06 0.54 ± 0.02 

Ri, (Mq) 0.56 ±0.06 0.50 ±0.04 0.53 ± 0.04 0.53 ± 0.04 0.53 ± 0.02 

[Fe/H] -0.056 ±0.2 0.04 ±0.13 ± 0.4 -0.05 ±0.17 -0.01 ± 0.09 

d (pc) ■ ■ • ■ • • ■ ■ • 303 ± 14 303 ± 14 

age (Gyr) ■ ■ • ■ • • ■ ■ • ~ 2.6 > 2 



References. — (1) IBrown et all (t201lD ; (2) IMuirhead eTaLl (I2012al> ; (3) 

IFabrvckv eTaLl l|2012al 'l 



were identified within AK' — 6.5 at a separation of 0"7 
or farther from Kepler-32. 

3.3. Kepler Light Curves 

Our light curve data come from the Kepler space tele- 
scope, which is conducting a continuous photometric 
monitoring campaign of a target field near the constella- 
tions Cygnus and Lyra. A 0.95-m aperture Schmidt tele- 
scope feeds a mosaic CCD photometer with a 10° x 10 ° 
field of view (IKoch et al.1 l20tnt IBorucki et all 1201 lbl) 
Data re duction and analysis is described by l Jenkins et ah! 
(|2010al lbl) and photometric and astrometric data were 
made publicly available as part of the 28 July 2012 public 
data release. We downloaded the data from the Mul- 
timission Archive at STScI (MAST), and we use the 
pi peline-corrected ligh t curves from Quarters through 
9 (jBatalha et"aT1l2012| ). 

These light curves were aggressively detrended to take 
out all astrophysical variation by first masking out all 
transit signals and then sequentially fitting a low-order 
polynomial to the time series data in 2 day chunks. These 
detrended data were then folded for each Kepler-32 d, e 
and f with the ot her planet signa l s mask ed according to 
the ephemeris of IFabrvckv et al. 1 (2012a). These folded 
transit signals were used for the FPP analysis. The 
binned photometric data are shown in Figure [2] with the 
trapezoidal fits used to assess the transit shape overlaid. 



3.4. False Positive Probability Analysis 

To calculate the relative probabilities of astrophysical 
false positive and bona fide transiting plan et scenarios, 
we follow the methodology of Morton (2012), which com- 
pares the observed transit shapes to those of simulated 
false positive populations, subject to the available obser- 
vational constraints. Although occurrence rates may be 
as high as 10% for planets in the radius bins correspo nd- 
ing to KOI-952.03, .04, and .05 ([Howard et al.lf2012h we 
assume a very conservative occurrence rate of 1%. 

Figure [3] summarizes these FPP results for all three 
signals. Despite our conservative estimates of planet oc- 
currence rate and without accounting for the fact that 
Kepler-32 b and c have already been confirmed to be 
coplanar transiting planets, which would further decrease 
the FPPs by a factor of about 10, we estimate the FFPs 
for KOI-952.03, .04, and .05 to all be less than 0.3%. We 
therefore consider all five photometric signals to be fully 
validated planet transits making Kepler-32 the richest 
known system of transiting planets around an M dwarf, 



TABLE 3 

Derived Planet Properties 



KOI 


Kepler 


P 1 


a 


R p 


rp a 
eq 






(days) 


(AU) 


(flffi) 


(K) 


952.05 


32f 


0.74296(7) 


0.0130(2) 


0.81(5) 


1100 


952.04 


32e 


2.8960(3) 


0.0323(5) 


1.5(1) 


680 


952.01 


32b 


5.9012(l) 


0.0519(8) 


2.2(1) 


530 


952.02 


32c 


8.7522(3) 


0.067(1) 


2.0(2) 


470 


952.03 


32d 


22.7802(5) 


0.128(2) 


2.7(1) 


340 



References. — (1) IFabrvckv et al.l l|2012al ') 
a assuming a Bond albedo a = 0.3 



and we assign the names Kepler-32 d, e, and f to the for- 
mer candidates KOI-952.03, .04, and .05, respectively. 
Verification of the transit signals allows us to take full 
advantage of the fortuitous alignment of the Kepler-32 
system with respect to our line of sight, and in the next 
section we look deeper into the details and implications 
of its specific configuration. 



4. THE Kepler-32 PLANETARY SYSTEM 

The increased accuracy with which we derive the stellar 
parameters for Kepler-32 (see TableEJ leads to more pre- 
cise planetary parameters, which we outline in Table [3l 
The planets of Kepler-32 have radii that are similar to 
values found in the Solar System, from below Earth size 
up to about 70% of Neptune. However, the planetary 
system of Kepler-32 has a remarkably distinct dynami- 
cal architecture in comparison. Figure |4] shows the rel- 
ative sizes of the planets and their orbits, along with 
labels denoting their periods, semimajor axes and period 
commensurabilities rounded to the nearest integer. The 
five planets of this system orbit within 0.13 AU from the 
star, or approximately one third of Mercury's semimajor 
axis. The outermost planet lies within a region where 
the stellar insolation is similar to Earth's, while at the 
other extreme, the innermost planet orbits only 4.3 stel- 
lar radii from the photosphere of Kepler-32. The three 
middle planets — termed e, b, and c in order of increas- 
ing semimajor axis — exhibit period ratios within 2% of 
a 1:2:3 commensurability. 

4.1. Tidal Evolution 

The proximity of the Kepler-32 planets to their host 
star imply that their dynamics have been significantly 
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Fig. 2. — Detrended, folded and binned Kepler lightcurves for 
the three planets that we validate with our FPP analysis. The 
trapezoidal fits from the analysis of [Morton (20l3) are overlaid. 
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Fig. 3. — The relative strength of each possible scenario in our 
false positive calculation — eclipsing binary (eb), hierarchical bi- 
nary (heb), background eclipsing binary (bgeb), background planet 
(bpl) and planet (pi) — is shown pictorially for our priors (upper left 
chart), likelihood (upper right chart), and final probability (lower 
chart) for each KOI-952.03 (top panel), .04 (middle panel), and .05 
(bottom panel). Observational constraints used in the calculation 
are stated in red (lower left) and the trapezoidal fit paramters are 
shown (upper right) in each panel. The quantity f p i^v ls t ne small- 
est value of the planet occurrence that would result in validation 
for the given transit signal. 
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Fig. 4. — Depiction of the Kepler-32 planetary system with the star and orbits drawn to scale. The relative sizes of the planets are shown 
at the bottom of the figure scaled up by a factor of 80 in relation to their orbits. 



altered due to tidal forces induced by the host star over 
the age of the system. Using an initial spin period of 10 
hours, a tidal dissipation factor of 100, a rigidity factor 
of 3 x 10 11 and a density profile similar to Earth such 
that the moment of inertia I ^ 0.5 MpR^, we derive the 
tidal locking timescales for the Kepler -32 planets to all 
be < lMyr (|Murrav fc DermottJ [2000h . This timescale 
depends linearly on the dissipation factor, Q, which is 
highly uncertain for exoplanets. However, for a reason- 
able range of Q values ranging up to 10 5 it can be ex- 
pected that the rotation periods of all the planets in the 
Kepler-32 system are equal to their orbital period barring 
spin-orbit resonances. 

Tidal forces from the host star will also damp the 
eccentricities of the plane ts in direct proportion to Q 
(|Murrav fe DermottJ l2000h . For Q < 10 4 the inner plan- 
ets will have no free (primordial) eccentricity. Though 
even for a low Q, the tidal damping time for Kepler- 
32 d, the most distant planet, is approximately 4Gyr, 
comparable to the age of the system. The free eccentric- 
ities of Kepler-32 b and c have recently been estimated 
to ha ve non-zero values from the phase of the T TV sig- 
nals (jLithwick et all 120121 IWu fc Lithwicl 120121) . This 
suggests high Q values for these planets, consistent with 
their moderate inferred densities (see § 14. 2[) . 

Contrary to the trend for K epler pairs to lie long - 
ward of mean motion resonance (Fab rvckv et al.ll2012bf i , 
Kepler-32b and c lie 1.1% shortward of resonance. 
The high Q values suggested above would limit dis- 
sipative mechanisms that may spread their orbits 
(|Lithwick fc Wd[2^ll[Ba"tvgin fc Morbidellill2012D . Also 
Kepler-32 b is near a 2:1 resonance with Kepler-32 e. 
However, our n-body simulation of § 14.41 suggests planet 



e has too small an effect to account for these effects. 

4.2. Mass and Density Estimates 

As mentioned above, Kepler-32 b and c were validated 
based o n the TTVs observed d ue to their mutual inter- 
actions (jFabrvckv et alj l2012af) . The amplitude of the 
TTV signals specify upper limits to their masses of be 
6.6 and 8.4 M®, implying densitie s of less than 3.4 an d 
5.7 gem -3 for b and c respectively (jLithwick et alJ 2012). 
However, the TTV phases are measured to be ~ 45° 
from and it implying a small free eccen t ricity . A cor- 
rection for this effect by IWu fc Lithwick (|2012f ) reduce 
the estimated masses for Kepler-32 b and c by factors 
of 0.51 and 0.45. With these corrections the masses 
and densities for Kepler-32 b and c are 3.4 and 3.8 M® 
and 1.7 and 2.6gcm -3 . The nominal masses for Kepler- 
32 b and c imply mass-radius relationships of M oc R lv , 
with 7 P = 1.5 — 1.9, similar to the value of 2.06 for the 
six Solar System pla nets bounded by Mars and Saturn 
(|Lissauer et al.ll2011l ). 

The above stated densities imply that Kepler-32 b and 
c are composed of a si gnificant amount of v olatiles. Using 
Equations 7 and 8 of iFortnev et al.l (|2007[ L we find that 
if Kepler-32 b and c had no atmospheres they would be 
expected to contain ~ 96% and ~ 56% volatiles. How- 
ever, given the equilibrium temperatures of Kepler-32 b 
and c a large fraction of their volatile content likely exists 
in the form of an atmosphere. 

4.3. Atmospheric Evolution 

The proximity of the Kepler-32 planets to their host 
star suggest significant atmospheric evolution due to 
evaporation, outgassing or both processes. The equi- 
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librium temperature of Kepler-32 f is ~ HOOK and its 
radius is measured to be 0.81 R®. For a planet this 
small with such a high equilibrium temperature the at- 
mosp heric mass fraction would have to be very small, ~ 
10~ 5 (Rog ers et al.ll20lU ). Using an extreme ultraviolet 
lumin osity of Kepler -32, L EU y ~ 10 26 6 dHodgkin fc Pvd 
119941 ) and following iLecavelier Pes Etangsl (|2007l ) using 
a conservative mass loss efficiency of euv — 0.1 we de- 
rive an atmospheric mass loss of ~ 10 8 gs -1 . Thus the 
timescale to lose its atmosphere is more than 100 times 
shorter than the age of the Kepler-32 system. We there- 
fore conclude the Kepler-32 f contains no atmosphere. 

Given the size and equilibrium temperature of Kepler- 
32 e its atmospheric mass fraction must also be small, 
M a /M p ~ 10~ 4 while the present day atmospheric mass 
loss rate is between 10 7 and 10 8 gs _1 . The timescale for 
the complete loss of the Kepler-32 e atmosphere is calcu- 
lated to be between 0.2 and 2 Gyr. Therefore Kepler-32 e 
must have lost a significant fraction of any atmosphere 
it started with. 

The total atmospheric mass loss for the other three 
planets is at least ~ 10~ 4 M s for reasonable choices of 
planetary mass. If these planets have relatively low- 
density cores (ice and rock) and started out with large 
atmospheres, they could have suffered considerable atmo- 
spheric evolution due to the heating by Kepler-32. Thus 
the observed sizes of the Kepler-32 planets are likely de- 
termined in part by the extreme ultraviolet and X-ray lu- 
minosity of their host star. However, the mass estimates 
from § 14.21 suggest that Kepler-32 b is 10% less massive 
than Kepler-32 c while being 10% larger and 25% closer 
to Kepler-32 hinting that the mass-radius relation for the 
Kepler-32 planets is not determined solely by a simple 
atmospheric evolution model. 

4.4. Kepler-32 Planetary System Architecture 

The physical characteristics of the Kepler-32 planets 
are summarized in Table [3] and the remarkably compact 
and orderly architecture of the system is shown schemat- 
ically in Figure |4j As mentioned above, three of the 
planets lie within 2% of a 1:2:3 period commensurabil- 
ity. Kepler-32 e and b have a period ratio of 2.038, which 
is 1.9% longward of commensurability, while Kepler-32 b 
and c have a period ratio of 1.483, or 1.1% shortward of 
commensurability. 

Planets within a mean motion resonance can stray a 
few percent from commens urability and maintain th e 
libration of resonant angles (|Murrav fc Dermottj |2000T) . 
However without detailed knowledge of the individual 
orbits it is not possible to determine with certainty if a 
planet pair is in a resonant configuration. Therefore we 
assess the significance of the near commensurability of 
Kepler-32 e, b, and c using a probabilistic argument. 

We randomly populate 5 planet systems with peri- 
ods between the inner and outermost planets in the 
Kepler-32 system enforcing separations larger than 2>/3 
for every pair of neighboring planets (IGladmanl 119931) 
and larger than 9 for chains of plane ts (Chamber s et al.l 
1996; Smi th et al.l 120091 iLissauer et al.ll2011l ) in units of 
mutual Hill radii. In this section a mass- radius rela- 
tionship of M oc R 2 06 (|Lissauer et al.|[20Tll) is adopted 
for consistency. The final ensemble of systems provide 
a baseline against which we can gauge the significance 



of the Kepler-32 architecture. The occurrence of period 
ratios in our randomly drawn sample that lie within ei- 
ther the observed period ratio of Kepler-32 e and b, or b 
and c, is fairly high, 11% and 12%, respectively. How- 
ever the fraction of systems arranged in a chain involving 
any combination of 1:2 or 2:3 near-commensurabilities 
is 1.9% ± 0.5%. This suggests that the architecture we 
observe today reflects the result of dynamical evolution 
rather than a result of chance. 

Given the compactness of the Kepler-32 system, 
we expand on the dy namical stability simulations of 
iFabrvckv et al . (2012a) by considering all five planets. 
Kepler-32 f is very likely on a circular orbit and lies > 30 
mutual Hill radii away from the nearest planet. It there- 
fore is not e xpected to have any bear ing on the stability of 
the system (|Smith fc Li ssaucr 2009) and is not included 
in the following analysis. 

Using the Mercury6 software (|Chambersj|1999D employ- 
ing the hybrid symplectic Bulirsch-Stoer integrator we 
integrate the Kepler-32 system for 30 Myr, or ~ 500 mil- 
lion orbits of the outermost planet. Our i nitial conditions 
were se t according to the ephemeris of IFabrvckv et ahl 
(|2012a|) with an initial time step of one tenth the orbital 
period of the innermost planet (in this case, Kepler-32 e). 
We began with zero eccentricity and random inclinations 
drawn uniformly between the minimum and maximum 
values that would produce a transit from our vantage 
point. We did not include any tidal effects. The eccen- 
tricities and inclinations varied stably with small ampli- 
tudes. Figure [5] shows the inclination evolution resam- 
pled at 100 kyr time steps and histograms of the eccen- 
tricities measured every 100 yrs over the 30 Myr simula- 
tion for the outer four planets. 

4.5. The Formation and Evolution of the Kepler-32 
Planetary System 

The long-term stability evident in our simulation of 
§ 14. 41 allows us to view the present-day architecture of the 
Kepler-32 planets as representing the state of the system 
at the end of its formation epoch, ~ 10 — 100 Myr after 
the formation of the star. We now investigate what can 
be learned about the formation of the Kepler-32 planets 
by looking backward from their end state. 

4.5.1. The Kepler-32 Protoplanetary Nebula 

Planets form from flattened disks of dust and gas cir- 
culating around protostars, and the masses and sizes of 
these protoplanetary disks as well as their lifetimes and 
accretion rates are con strained from observations (see 
I Williams fc Ciezall2011l) . The surface density profiles are 
typically assumed to follow a power law form cx a 1 . 
(Note: Here and throughout we use a to denote both the 
orbital radius of circumstellar material and the semima- 
jor axis of the plan ets. The Kepler plane ts typically have 
low eccentricities (|Wu fc Lit hwick 2012) minimizing am- 
biguity.) 

The value of 7 is often taken to be 3/2 based on 
estimates of the minimum mass so lar nebula (MMSN; 
lWeidenschillinglll977t lHavashill 19811 ) . The MMSN is con- 
structed by smoothing the planets of the Solar System 
over their respective domains and correcting for Solar 
abundance to recreate a minimal surface density profile 
from which the Solar System could have formed. Mod- 
ern observations constrain the surface density profile in 
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Fig. 5. — Time evolution of the inclinations sampled every lOOkyr 
and the distributions of eccentricities for the Kepler-32 planets 
sampled every 100 yrs from our n-body simulation. Kepler-32 f was 
not included in the simulation due to its dynamical isolation from 
the rest of the planetary system (see text). 



the outer disk (> 40 AU) to be on average shallower than 
this. 7 » 0.4-1.0 (lAndrews fc Willia n9l2007t lls ella et al.l 
120091: lAndrews et al.ll2009D . while the form of the surface 
density profile in the inner disk remains largely uncon- 
strained. 

We can construct a "minimum mass protoplanetary 
nebula" for Kepler-32 in a similar fashion. Since we 
do not have mass estimates for all the planets in the 
Kepler-32 system we use the mass- radius relationship fo r 
Solar System planets M oc i? 206 (|Lissauer et al.ll2011j) . 



consistent with the nominal masses from § 14.21 to within 
50%. Assuming a gas-to-dust ratio of 100 and a rocky 
core mass fraction of 50% implies that the protoplanetary 
nebula of Kepler-32 contained ss 3 Mj up of gas between 
0.013 and 0.13 AU of the host star. 

This translates to a very high surface density, sa 5 x 
10 5 gem -2 , in the inner regions of the disk that is incon- 
gruous with the masses and surface density profiles in the 
outer regions of disks as measured with (sub)-millimeter 
interferometrv (|Andrews fc Williams! 120071 : llsella et al.l 
120091 : lAndrews et al.l I2009T ) . Indeed, if this surface den- 
sity is used to scale the median best fit values of the disk 
radius and surface density profile index for stars less than 
1 M Q in the lAndrews et al.l (|2009( ) sample— R d = 126 AU 
and 7 = 0.9 — the total mass of the Kepler-32 protoplan- 
etary nebula would be orders of magnitude greater than 
the star! This implies one of three possibilities: the disk 
surface density deviated significantly from a single pow- 
erlaw in the inner regions with a larg e pile- up of material 
near the disk inner radius (see, e.g., iChiang fe Laughhnl 
120121) : the material that formed the planets came from 
elsewhere in the disk; or the planets themselves formed 
elsewhere and were transported to their present-day lo- 
cations. 

4.5.2. Oligarchic Growth and the Isolation Mass 

The gravitational influence of the host star poses an- 
other barrier to super-Earth-sized planets forming within 
a tenth of an AU. The amount of material accessible 
to a growing planetary embryo during the oligarchic 
phase of growt h is estimated by the "isolation mass" 
(|Lissauerill987D . This quantity is derived by calculat- 
ing the amount of mass available to an oligarch within a 
disk of planetesimals. We choose the radius of gravita- 
tional influence to be 4 Hill radii based on the approx- 
imate spacing of isolated o ligarchs in numerical simula- 
tions (jKokubo fc Idalll998D to obtain 
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where a is the radial distance from the star, M* is the 
mass of the star, and S(a) is disk surface density pro- 
file. To calculate values for M iso (a), we again use the 
median surface dens ity profile of the M stars in the 
lAndrews et~aT1 (120091) sample, the median values of the 
M stars in the llsella et al.l ( 120091) sample w ith 7 = 0.5 
and R d = 260 AU, and also the lEisnerl (|2012D model with 
a median disk radius of Rd = 250 AU and a surface den- 
sity profile with fixed 7 = 1.37 based on theoretical ar- 
guments. 

Figure [5] shows that the isolation mass computed for 
these surface density models are > 3 orders of magni- 
tude too small to account for the mass in the Kepler- 
32 planets. While the planetary embryos can undergo 
an additional phase of growth through subsequent merg- 
ers t hat increase their mass by a n order of magnitude 
(e.g.. IChambers fc Wetherililll998f ). it would be difficult 
to achieve an addit ional factor of ~ 100 growth during 
this phase (see also lLissauerll2007l) . The total amount of 
growth in this late phase growth is limited to the total 
amount of mass present in the disk at the outset, which 
is addressed in the previous section. 
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of the dust sublimation radius for pre-main sequence ages of 2 Myr and 10 Myr IIBa raffc ct al. 1998), the approximate location of the ice 
line during the planet forming epoch (blue shaded region, Kennedy & Kcnvon 2008), and the stellar radius of Kepler-32. The discovery 
space of the Kepler mission through quarter 6 (from which this sample is drawn) is roughly delineated with the thin dashed lines. 



The migration of planetesimals from outside this re- 
gion also offers a p ossibility for late time gr owth of the 
Kepler-32 planets ( Hansen fc Murravl 120121 ). However 
the gravitational stirring of the planetesimals can cause 
radial displacements of the planetary embryos away from 
period commensurability. The magnitude of this effect 
has recently been used to place limits on the amount of 
mass contained in orbit-crossing planetesimals at the end 
stages of planet formation for closel y spaced, low incl i- 
nation multi-planet Kepler systems (Moore c t al.l [20121 ) . 
These limits of ~ 1 Mq are too low to contribute signifi- 
cantly to the discrepancy between the isolation mass and 
the final mass of the Kepler-32 planets. 

This corroborates the above conclusions that the 
Kepler-32 protoplanetary disk either contained a dispro- 
portionately large amount of mass within 0.1 AU, or that 
the formation of the Kepler-32 planets followed a more 
complex formation and evolution history. 

4.5.3. The Dust Sublimation Radius 

Planets form from dust in protoplanetary disks. But 
at small orbital radii the stellar irradiation raises the 
dust temperature beyond the sublimation point, clear- 
ing these regions of the raw materials for planet forma- 
tion. Low-mass stars such as Kepler-32 begin pre-main 



sequence evolution with luminosities much larger than 
their main sequence values and take > 10 Myr to con- 
tract to the point when hydrogen burning begins. The 
sublimation radius at these early epochs is thus further 
from the star than one would conclude from the present- 
day luminosity of the star. 
The dust sublimation radius is given as 



Srai = 0.034 ("9°Y 

V sub I 




(3) 



(jlsella et al.l I2009D where T su t, is the sublimation tem- 
perature of the dust, L* is the stellar luminosity, and e 
is a measure of the cooling efficiency of the grains with 
larger grains havi ng higher efficienc ies leading to smaller 
sublimation radii (jlsella et al.ll2006D . Figure [7] shows the 
evolution of R su b as a function of pre-main sequence age 
of a 0.5 Mq star bas ed on the stellar evolution models of 
IBaraffe et al] (Il998l) . 

Even if it is conservatively assumed that the dust is 
in large grains (e = 1) and composed of p ure iron with 
a sub limation temperature T su b = 1908 K ( Pollac k et al.l 
Il994f ). the sublimation radius does not move inward of 
Kepler-32 f's orbit until after 9.8 Myr of pre-main se- 
quence evolution. For smaller grains made of silicates 
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Fig. 7. — Dust sublimation radius of Kepler-32 as a function of 
pre-m ain sequence age according to the models of Bara ffe et al.l 
(1998). The solid curve is for large iron grains (T su f, = 1908 K 
and e = 1) the dashed curve is for large silicate (olivine) grains 
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(olivine) grains within a lower density environment (T su (, = 1570 K 
and e = 0.58). All curves cross the present day semimajor axis of 
Kepler-32 f after the typical lifetime of a protoplanetary disk. 



and for less dense environments the dust sublimation ra- 
dius lies beyond the semimajor axis of Kepler-32 f un- 
til after ~ 50 Myr, and beyond the present day orbit of 
Kepler-32 e for ~ 3 Myr of pre-main sequence evolution. 
This poses a serious challenge to in situ formation par- 
ticularly since the dusty disks from which planets form 
only last < 10 My r before they are drained , disrupted, or 
evaporated away ([Williams fe Ciezall201lD . 

4.5.4. In Situ Formation vs. Formation then Migration 

The above investigation into the structure of the 
Kepler-32 planetary system as well as the history and 
evolution of the host star present meaningful constraints 
on its formation. To form Kepler-32 f and e where we 
view them today requires that they formed from plan- 
etesimals scattered inward from the wider orbit plan- 
ets. This kind of formation scena rio would be dif- 
ferent than the theory described by Hansen fc Murray! 
(|2012t) where scattered planetesimals accrete onto pre- 
existing oligarchs from an early phase of planet forma- 
tion. Rather Kepler-32 f and e would have had to have 
formed from scattered planetesimals in a region largely 
devoid of solid material. It is not clear how feasible this 
would be, and requires an investigation beyond the scope 
of this work. An additional challenge for in situ forma- 
tion of the Kepler-32 planets is to explain the near period 
commensurabilities of Kepler-32 e, b and c with in situ 
formation that we have shown to be statistically signifi- 
cant beyond pure chance. 

On the other hand, a scenario in which the Kepler-32 
planets formed further out in the disk and then migrated 
to their present locations offers a natural explanation for 
1) the large amount of planetary mass within 0.1 AU of 
this 0.54 Mq star, 2) the position of Kepler-32 f within 
the dust sublimation radius of the star for ~ 10 Myr 
of pre-main sequence evolution, and 3) the near period 
commensurabilities between three of the planets. The 
high volatile content inferred for Kepler-32 b and c is 



also indirect evidence that these planets may have ac- 
quired volatiles from beyo nd the snow line at ~ 1 AU 
(|Kennedv fc Kenvonll2008f) . The difficulties in applying 
this concept lie in understanding the mechanisms of mi- 
gration through a disk of gas and planetesimals. While a 
full, rigorous analysis of these possibilities is also beyond 
the scope of this article, we look to some order of magni- 
tude calculations of different migration mechanisms for 
insight. 

4.5.5. Migration 

The compact and coplanar architecture of the Kepler- 
32 system favors migration through a disk rather than 
through planet-planet interactions. Interactions with ei- 
ther a gaseous disk or a disk of planetesi mals can cause 
radia l migration of planetary orbits (see Kie v fe Nelson! 
120121 for a review). Planets massive enough to open a 
gap in the protoplanetary disk are coupled to the viscous 
evolution of the disk while less massive planets migrate 
more quickly due to larger disk torques. 

Using a ge neralized criterion for gap formation in a 
viscous disk ([Crida et al.l [2006. Eq. 15) we find that 
none of the Kepler-32 planets would have opened a gap 
in their protoplanetary disk assuming a constant disk 
flaring term a/H = 0.05 where H is the scale height, 
a viscosity parameter a = 10~ 3 and using the nomi- 
nal masses for Kepler-32 b and c (§ I4.2|) and M p oc i?p 06 
otherwise. We therefore restrict our discussion of migra- 
tion through a gaseous disk to the "type I" mechani sm 
(|Goldreich fc Tremaind[l979t iLin fe Papaloizoulll979! ) 

Assuming a protoplanetary disk with with a power- 
law surface density profile, S = Erjo -7 > and a midplane 
temperature profile, T oc a~^ T , we esti mate the migra- 
tion ra te through a gaseous disk following Kie v fc Nelson! 
(2012) 



da = _2_ 

dt m ri V GM+ 



(4) 



where p subscripts refer to planet quantities, and r to t is 
the total torque on the planet from Linblad, corotation, 
and horseshoe torques 



rtot = [-1.36 + 0.627 + O.430H 
(also see lMasset et al.l [2006). and 



1.36(3/2- 7 )] r„ (5) 



To 
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dp J 
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As a baseline model we use a disk surface density pro- 
file with 7=1, which extends from the present day lo- 
cation of Kepler-32 f out to 200 AU, contains a total of 
10% of the stellar mass and has a midplane tempera- 
ture profile with fix = 3/4. This disk model yields type 
I migration rates da/dt ~ —0.07,-0.3,-0.7,-0.9 and 
-1 AU (lO-^yr)" 1 for Kepler-32 f through Kepler-32c, re- 
spectively with the migration rate scaling as da/dt oc 

m p Eoap^ 2-7 . For these calculated rates, the timescale for 
these planet s to have migrated from bey ond the snow line 
at ~ 1 AU (|Kennedv fc Kenvonl I2008D to their present 
day locations is short, r mig ~ 10 4 years, in comparison 
to typical disk lifetimes. 

The increasing mass of the Kepler-32 planets as a func- 
tion of semimajor axis naturally produces convergent mi- 
gration for type I torques in a smooth disk, a necessary 
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but not sufficient condition for resonant capture. To es- 
timate the probability for resonant capture we perform 
an order of magnitude calculation comparing the reso- 
nance crossing times to the libration time for the reso- 
nances. The resonance widths and libration times are 
estimate d using the pendulum mo del for an interior res- 
onance (jMurrav fc Dermottl 120001 §8.6-7). For the 2:3 
resonance of Kepler-32 b and c, the libration period is 
Tm, ~ 10 years and and the width is estimated to be 
Aa res ~ 0.009 AU using the mean eccentricity of Kepler- 
32 b from the simulations of § 14.41 With a convergence 
rate of 0.2 AU (10 4 yr) _1 we find that the resonance cross- 
ing time is roughly 40 libration periods. Following the 
same procedure for the 1:2 resonance of Kepler-32 e and 
b, we find that the resonance crossing is of order 10 li- 
bration periods. 

Shallower surface density profiles produce slower con- 
vergence rates, increasing the probability of resonant 
capture while steeper profiles produce faster convergence, 
shortening resonance crossing times and decreasing the 
probability for capture. For 7 = 1/2 we find that the res- 
onance crossing times are on the order of 10 3 libration 
periods for both resonances and for 7 = 3/2 these times 
are shortened to < 1. The resonance crossing times in 
units of libration periods is a weakly increasing function 

of semimajor axis for our simple model, going as a}J 4 . 

Therefore we find that the Kepler-32 planets could 
have been captured into resonance while migrating 
through a gaseous disk if the surface density profile had 
7 < 1, consiste nt with modern obse rvations of proto- 
planetary disks ([Andrews et alJl2~009h . The probability 
for capture is a weak function of semimajor axis, and 
it is thus also feasible for the Kepler-32 planets to have 
been caught in resonance further out in the disk and 
have moved inward in lockstep (e.g. ICresswell fe Nels"onl 
2008). How the planets may have stopped their migra- 
tion before falling onto the star remains an outstanding 
question in this scenario. 

It is also possible that the Ke pler-32 planets mi- 
grated through a planetesimal disk (jLevison et al .1120071) 
rather than the gas dominated disk required for type 
I migration. We calculate the order of magnitude 
ra tes for planetesimal migra tion using the equations 
of iBromlev fc Kenvonl ([201 If ) assuming a gas to dust 
ratio of 100 and that all solids are in the form of 
planctesimals. We find migration rates of da/dt ~ 
-0.44, -1.1, -1.7, -2. 1,-3. 2 AU (10 s yr)" 1 , roughly 4 
orders of magnitude slower than type I. This slow mi- 
gration rate would favor resonance capture and may cir- 
cumvent the need for an efficient stopping mechanism. 
However, it may be difficult to migrate these planets in 
lock step from much further out in the disk as the in- 
ner planets could stir the planetesimal disk reducing the 
migration e fficiency at that location for the next furthe r 
out planet (|Kirsh et alj|2009t IBromlev fc Kenvonll20ill ). 
Also if Kepler-32 b and c were migrating at such slow 
rates, it is likely that they would have ended up in a 
1:2 resonance rather than the more compact 2:3 reso- 
nance. These details may limit the distances over which 
planetesimal migration could have acted to produce the 
observed architecture of Kepler-32. 

From these analyses, migration through a gaseous disk 
produces a more favorable situation for transporting the 



Kepler-32 planets from far enough out in the disk to 
reconcile the difficulties with in situ formation outlined 
above. However, further investigation into the migration 
of multi-planet systems through both gas and planetesi- 
mal disks is needed to draw more definitive conclusions. 

5. DISCUSSION 

The analyses of the preceding sections provide evidence 
that the Kepler-32 planets formed further out in their 
protoplanetary disk from where we see them today and 
migrated convergently into their present locations. The 
high inferred volatile content of Kepler-32 b and c, and 
the order of magnitude calculations for the migration 
rates of the Kepler-32 planets suggest that these plan- 
ets have formed beyond the snow line in the presence 
of gas. If true, the Kepler-32 planets would have neces- 
sarily formed within ~ 10 Myr, the ti mescale over which 
gas s urvives in protoplanetary disks ([Williams fc Ciezal 
1201 If ) - This formation history stands in contrast to the 
formation of the terrestrial planets in the Solar System, 
which are commonly thought to have formed i n situ in a 
gas- fr ee environment on a 100 Myr timescale ([Wetherilll 
1990). The possibility to constrain the timescale for the 
formation of the Kepler-32 planets motivates further in- 
vestigation into migration of multi-planet systems. 

Our conclusions about the formation of the Kepler- 
32 planets rely on a detailed characterization of a single 
system not possible for the full ensemble of Kepler's M 
dwarf planets. However, both the Kepler-32 star and 
its planets are representative of this full ensemble giving 
important clues regarding the formation of full sample of 
Kepler M dwarf planet candidates. 

5.1. The Ensemble of Kepler M Dwarf Planets 

Our sample of M dwarfs was drawn f rom the cata- 
log of Kepler Objects of Interest given by Batalha et al. 
(|20T1) using a color magnitude cut oi Kp > 14 and 
K p - J > 2 (|Mann et all 120121 ). This results in 5499 
stars from the KIC catalog observed with Kepler, how- 
ever there are only 4682 with finite, non-zero photo- 
metric precision values for at least one quarter between 
Q0 and Q6. Since the Kepler sample is magnitude lim- 
ited, the sample is skewed toward the most massive M 
dwarfs, peaking around 0.5 Mq (see Figure ITOl) . Two gi- 
ant planet candidates are present in this sample. One 
is a bona fide gi ant planet confirmed with radial velocity 
data fKOI-254: [Johnson et al.lf2C>T2h . The other, KOI- 
1902, we found to be a false positive based on its transit 
profile. We scrutinized the light curves of several of the 
other apparent outliers by hand. We reject two KOIs in 
the list, KOI-531 and KOI-1152, on account of deep sec- 
ondary transits, and two others, KOI-256 (Muirhead et 
al. in prep.), and KOI-1459 based on their transit pro- 
files. The final distribution of M dwarf planet candidates 
are shown as gray circles in Figure [5] These remaining 
candidates are expected to have a high fidelity allowing 
us to treat them as a sound, statistical ensemble. 

The distribution of M dwarf planet candidates in Fig- 
ure [6] contains 100 planets in 66 systems. There are 48 
single planet systems, 7 double planet systems, 7 three 
planet systems, 3 four planet systems, and 1 five planet 
system — Kepler-32. Therefore 18 systems or 27% of the 
M dwarf KOIs are multi-transit systems. The giant, 
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KOI-254 is a single candidate system that constitutes 1% 
of all planets in the Kepler M dwarf sample and about 
2% of all planet hosting M dwarfs. This planet is anoma- 
lous and we do not consider it as part of the ensemble 
hereafter. 

The main distribution of M dwarf planets appears to 
follow a trend of increasing mass (radius) as a function 
of semimajor axis. Without accounting for biases, we 
derive the relationship, M oc a o e . However, for a given 
signal to noise threshold the lower envelope of planet 
candidates is expected to follow 

i? i „„ = (SNRa p ) 1 /2 i? 2/V/ 4 (^) 1/4 (7) 

where SNR is the signal to noise threshold, a p is the 
photometric precision, R+ is the stellar radius and n is 
the number of transits measured. For our adopted mass- 
radius relationship this translates to Mu m cx a 5 . The 
detection limit for an observational baseline of 6 quar- 
ters and a SNR threshold of about 5 for a typical M1V 
star is shown in Figure [6] as the dashed line. This sug- 
gests that the apparent trend at the lower envelope of the 
distribution is due to observational biases. Also shown 
in Figure [5] is the semimajor axis corresponding to Ke- 
pler's observing time baseline over the first 6 quarters 
of data over which t he planet candidates were selected 
(jBatalha et aL|[2012h . 

Although we cannot perform as detailed analysis on the 
Kepler M dwarf planet ensemble as we have for Kepler- 
32, we note that positions and inferred masses of the 
planet candidates as a whole also imply very large dis k 
masses within 0.1 AU (also see lChiang fc Laughlinll2012f ). 
Figure [8] shows the reconstructed mean protoplanetary 
disk masses for the 18 multi-transit M dwarf KOIs as a 
function of 7, the 1 a spread in the derived masses for any 
given value of 7 is approximately 0.3 dex (shown as data 
point and error bar). Even for 7 = 3/2, an appreciable 
fraction of the disks would be gravitationally unstable 
and therefore unlikely to produce the observed planets 
where we see them. 

We also analyze the present-day locations of the Kepler 
M dwarf planets in terms of the sublimation radius of 
their host star at 10 Myr of pre-main sequence evolution. 
We find that between 5 and 14% of the total number of 
planet candidates fall within the dust sublimation radius 
of their host stars for large iron grains and 1 /im silicate 
grains, respectively, and could not have formed in place. 

5.2. Kepler-32 as a Representative of the Kepler M 

Dwarfs 

The Kepler-32 planets span the main distribution of 
the Kepler M dwarf planets seen in Figure [5] and we plot 
in Figure [§] the locations of these planets in relation to 
the full distribution as a function of planet radius and 
semimajor axis. In both parameter spaces the Kepler-32 
planets fall well within the main distribution of planet 
candidates. 

To further explore how representative the Kepler- 
32 planetary system is of the full sample of Kepler 
M dwarf planet candidates, we create an ensemble of 
planetary systems with the Kepler-32 system specifica- 
tions oriented randomly on the sky with mutual incli- 
nations drawn randomly from a Rayleigh distribution 




-2t_^ , 1 , , , , 1 , , , , 1 , , , 

0.5 1.0 1.5 2.0 

Y 



Fig. 8. — The inferred disk mass in units of the stellar mass 
for all the M dwarf multi-transit systems in the Kepler sample as 
a function of assumed surface density profile index 7, and disk 
radius, i?disk- The data point and error bars show the mean and 
standard deviation of one such calculation for 7 = 1.5 and -Rdisk = 
150 AU to give a sense of the distribution of disk masses derived in 
each calculation. Disks with masses > 0.2 M* are expected to be 
gravitationally unstable. 



(jLissauer et al.l 1201 If ). We find the transit multiplic- 
ity fractions are best reproduced with a spread in mu- 
tual inclinations of 1.2° ± 0.2°, consistent with values 
for t he entire Kepler plane t candidate ensemble (1.0°- 
2.3°: iFabrycky et al.ll201 2b). The discrepancies between 
the real and simulated distributions of transiting systems 
are largest for the single transit systems (55% simulated 
vs. 73% real). This would be expected for a situation 
where the typical planetary system contains less than 5 
planets or is less compact. However, we note that the 
observed transit multiplicity of Kepler M dwarf systems 
can be recreated remarkably well assuming the case that 
all planetary systems are exact clones of Kepler-32. 

The distribution of stellar masses, radii, and metallic- 
ities for the Kepler M dwarf sample are shown in Fig- 
ure 1101 A narrow range of stellar mass and radius is 
seen peaked around 0.5 in solar units due to the mag- 
nitude limit of the sample and is shown in contrast to 
the present da y mass function of single stars (PDMF; 
Chabrier 2003) in the top panel. Kepler-32 is seen to be 
representative of the full sample in all quantities. 

5.3. Planet Occurrence 

The planet occurrence rate for M dwarf planets has 
been calculated to be about 0.3 for planets with rad ii 
> 2 i? e and periods less than 50 days (jMann et all 2012). 
However, only 28 planets of the 100 total in our Kepler M 
dwarf ensemble satisfy these criteria. Therefore the total 
planet occurrence rate for short period planets around M 
dwarfs is much higher than this number. The detection 
of planet signals within the Kepler M dwarf sample is 
not uniform causing a systematic uncertainty to planet 
occurrence estimations. Since we have no way to correct 
for this currently, we ignore this effect which will result 
in a lower limit of the occurrence rate. 

The number of planets per star in our sample is esti- 
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Fig. 9. — The distributions of Kepler M dwarf planet candidates 
as a function of planet radius and semimajor axis (blank histogram; 
not including KOI-254, 256, 531, 1152, 1459, or 1902) in compari- 
son to the Kepler-32 planets (hashed histogram). 



mated as 



En* i 
j ' Pi,j 



(8) 



where j is the index of the sum over all stars around 
which planet i c ould be detected, (see, e.g., 

Ho ward et aT1l2012D . and 



Pi,3 



(9) 



is the geometric probability of detecting planet i around 
star j if eccentricity is negligible and Rp/R* <C 1. Since 
the KIC stellar radii values for stars with T c g < 4500 K 
are unreliable, we randomly sample the distribution of 
stellar radii in Figure [10] estimated from near infrared 
spectra for this calculation. Evaluating Equation [8] over 
our sample yields / = 1.0 ± 0.1 planets per star where 
we use a binomial error estimate. 

Though we do not perform statistical validations for 
all the Kepler M dwarf planets, our sample is expected 
to have a higher fidelity than the total ensemble of Ke- 
pler transit signals since we culled our sample by hand. 
The planet occurrence derived from our sample assum- 
ing 90% fidelity, /go, is calculated by evaluating equa- 
tion [5] repeatedly for 90 randomly drawn planets of the 
100 total. We find f 90 = 0.9 ± 0.1 where the stated er- 
ror is the standard deviation of the distribution of 1000 
realizations. As a futher check, we limit our sample to 



20 
15 



10 



- I I I 

: <M.> = 0.51 M sun 
" ° M = 0.07 


i,im 














0.1 



0.2 



0.3 0.4 0.5 
Stellar Mass (M sun ) 



0.6 



0.7 



20 
15 



10 



I I I I ' ' 




: <R.> = 0.49 R sun 




7 o R = 0.07 I — 










:m : 



0.1 0.2 0.3 0.4 0.5 
Stellar Radius (R sun ) 



0.6 



0.7 




0.0 
[Fe/H] 



0.4 0.6 



Fig. 10. — The distribution of stellar masses and radii for our 
sample of M dwarf planet host stars are sharply peaked around 
0.53 Mq and 0.51 Rq, respectively. The present day mass function 
of single field stars in the Galaxy ijChabricr 2003) normalized to 
the number of stars in our sample is shown as the dotted line in 
the top panel for reference. Metallicities for our sample fall mostly 
between —0.4 and 0.4. Kepler-32 values denoted by the vertical 
dashed line all fall near the center of these distributions. 

the 28 planets with R p > 2 and P < 50 d to obtain 
f = 0.26 ± 0.05 a t the lower end of the estimations of 
iMann et al.l (|2012| ). This may be further evidence that 
our estimations are conservative. 

Counting only one planet per system in equation [8] 
gives the occurrence of stars with at least one planet, 
0.51 ±0.07 (not including KOI-254), which together with 
our estimate of the total occurrence gives the average 
number of planets per system as ~ 2. Stars with masses 
characteristic of the Kepler M dwarf sample are a fac- 
tor of ~ 1.8 t imes more common than stars of 1 Mq 
( Chabrier 2003) , and the planets of the Kepler-32 system 
are representative of the planets that form around these 
stars. Therefore the insights gleaned from the Kepler-32 
system both show us where to look for additional planets 
in the Solar Neighborhood, and provide a template for 
understanding the formation of the ubiquitous compact 
planetary systems throughout the Galaxy. 

5.4. Summary and Future Directions 

We present a detailed analysis of the Kepler-32 plane- 
tary system which offers the rare circumstance of 5 tran- 
sit signals. While two of the planets have been previously 
validated through evidence of their mutual gravitational 
interactions, we validate the remaining three transit sig- 
nals probabilistically using observations from the W. M. 
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Keck Observatory as constraints. This validation makes 
Kepler-32 the richest system of transiting planets known 
around an M dwarf. 

Kepler-32 has a markedly compact architecture. All 
five planets orbit within one third of Mercury's distance 
from the Sun, with the closest planet orbiting only 4.3 
stellar radii from the Kepler-32 photosphere. The three 
middle planets lie close to a 1:2:3 period commensurabil- 
ity that is unlikely to be the result of chance. 

Our refined stellar parameters improve the derived 
planetary characteristics, and aid in reconstructing this 
system's formation history. Several pieces of evidence 
from our analyses indicate that the Kepler-32 planets 
did not form where we see them today: 

• the dust sublimation radius of Kepler-32 lying out- 
side the present day semimajor axis of Kepler-32 f 
for longer than a typical protostellar disk lifetime 

• the extremely high surface densities inferred by as- 
suming in situ formation 

• the limited range of gravitational influence for plan- 
etary embryos located so close to their host star 

• the unlikely arrangement of three planets in the 
system near a 1:2:3 period commensurability 

• the high volatile content of Kepler-32 b and c. 

This conclusion necessitates planet migration through 
a disk, and our order of magnitude calculations for the 
migration rates of the Kepler-32 planets embedded in a 
typical protostellar disk suggest the presence of gas. If 
true, this would limit the formation time of the Kepler-32 
planets to < lOMyr — the known timescales over which 
gaseous disks survive. 

Kepler-32 is found to be representative of the full sam- 
ple of 66 Kepler M dwarf host stars, and the Kepler-32 
planets span the mid-line of the distribution of 100 Ke- 
pler M dwarf planet candidates in radius-semimajor axis 
space. Although we are unable to treat each system in 
this ensemble with the same care as for Kepler-32, we 
show that similar analyses applied to the ensemble give 
consistent results to those derived for Kepler-32. Thus 
the formation scenario deduced from Kepler-32 offers a 
plausible blueprint for the formation of the full sample 
of Kepler M dwarf planets. 

We select out 4682 M dwarfs from the Kepler Input 
Catalog that have been observed with Kepler and use 
their observational parameters to derive the planet oc- 
currence rate of Kepler M dwarf planet candidates. We 
confirm that within the completeness limits of the first 6 
quarters of Kepler data, the M dwarf planet candidates 
have an occurrence rate about 3 times that of solar-type 
stars, while the occurrence rate of all candidates around 
M dwarfs is 1.0±0.1. We expect the fidelity of our culled 
sample to be above 90%. Thus the compact systems of 
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planets around the Kepler M dwarf sample are a major 
population of planets throughout the Galaxy amplifying 
the significance of the insights gleaned from Kepler-32. 

At the time of this writing, there are only 37 planets 
confirmed to ex ist around 24 M d warfs in the Galaxy 
(exoplanets.org. [Wright et al.ll201l1) . The K epler space 
telescope has revealed 100 planet candidates around 66 
M dwarfs from which we draw our statistics. It would be 
of great benefit to the study of M dwarf planet formation 
to explore and validate a larger statistical sample of M 
dwarfs such that comparisons can be dra wn against th e 
detailed analyses of solar-type stars (e.g., lYoudinlfeOllh . 
The continued monitoring of the current sample is also 
important as this will reveal trends in planet occurrence 
as a function of orbital period for the smallest planets. 

Mass measurements of the growing numbers of con- 
firmed M dwarf planets will also play an important role 
in interpreting their origins. This can be achieved with 
nearby M dwarfs using precision radial velocity measure- 
ments, or by using alternative techniques such the am- 
plitude of transit timing variations. 

Lastly, direct imaging of the inner few AU of nearby 
protostellar disks will be possible in multiple wavebands 
with the Atacama Large Telescope Array (ALMA). The 
modeling of this emission may be the most direct way 
to constrain where in the protoplanetary disk compact 
planetary systems form. 
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